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1 Introduction 



The real photon production and dilepton production in a hot plasma are intensively discussed 
subjects of the current literature. One aspect of this discussion is the behaviour of the produc- 
tion rate near zero invariant mass of the lepton pair. Numerical simulations find a vanishing 
rate, while higher loop perturbative p[ and HTL-improved calculations find diverging result 
in this limit. The possibility that the two-loop order enhancement may be a consequence of 
infrared (IR) singularities was studied in recent papers [Q, where also an effective description 
was suggested to overcome the problem. In a resummation method was worked out for 
ladder diagrams, and it was applied to the photon and gluon production rate. 

Resummation in perturbation theory becomes unavoidable, in general, when we encounter 
IR divergences that enhance the contribution of higher order loops, thus overriding the cou- 
pling constant suppression. A well known example is the on-mass-shell IR divergence of the 
propagator due to a radiative mass shift. This divergence is cured by resumming IPI diagrams 
with the help of Schwinger-Dyson (SD) equations. This example also illustrates that resumma- 
tion may be effective only for some specific momentum range, while in other regions its effect 
might remain negligible. From this point of view, also finite-temperature Debye screening is a 
consequence of IR divergences. 

In non-equilibrium systems there appear new sources of IR divergences. When working in 
real time one recognizes that higher loops may give contributions that grow as some power of 
the elapsed time (~ t"" where n is the loop order). This would clearly restrict the applicability 
of perturbation theory to rather short times, so we need a resummation. It turns out that 
resummation leads to Boltzmann equations describing the temporal evolution of the propagator. 
In Fourier space the corresponding phenomenon is the appearance of pinch singularities 0, |^ 
in two-particle intermediate states with vanishing momentum difference. The resummation of 
these divergences, which involves ladder diagrams, leads, again, to Boltzmann equations 0. 

In this paper we shall argue that there are further momentum ranges where we can observe 
pinch singularities. When these momentum ranges are relevant, we have to perform appropriate 
resummation also there. To see the origin of these divergences we have to analyse the corre- 
sponding pinch singularities, as will be done in Section ^. The core effect of pinch singularities 
is that in a pair of propagators one of them can force the other to its mass shell; for instance 
in G r{j))G A^q) = GR{p)ig{q) + G r{p)G R^q) the first term yields the singular effect, the second 
turns out to be regular. If we parametrise the momenta carried by the particles asp = Q + K/2 
and q = Q — K/2, then we can find for the contributions two formally equivalent expressions: 

6{{Q-K/2f-m^) ^ 6{{Q-K/2f-m^) ^ I 6{{Q - K I2f - m^) 
R{P)m) [Q + K/2Y-m^ 2QK 2 Q^ + K^/A-m? ' ^' 

(The Landau prescription is omitted to simplify this introductory argument.) Thus the product 
of the propagators becomes divergent either when 2QK = (this case of the resummation leads 
to the Boltzmann equation), or when + K^/A — rn? = 0. This second is the new kinematical 
situation studied in the present paper. 

After identifying the IR divergences, we have to design an appropriate resummation scheme 
for them. In order to perform ladder resummation we can use the linearized Kadanof-Baym 
(KB) equations 0. These come from a procedure (for details see Sections ^]3| and |), which 
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consists of writing two SD equations for each of the Wigner-Fourier-transformed propagators, 
next taking their sum and difference, and finally linearizing them. They read, schematically: 



2QKG{Q,K) = {A.*G){Q,K), 2{Q^ - m^)G{Q,K) = {A+* G){Q,K), (2) 

where A± denote the linear operators acting on the propagators. The equations are diagonal 
in K, which is the momentum associated with the average coordinate of the propagator. The 
coefficient of G vanishes on the left-hand side in regions of pinch singularity, while the expression 
is finite on the right-hand side. The perturbative solution diverges at every order in any 
such regime, but the complete solution of the differential equations is regular. This is the 
resummation realized with the help of SD equations. Using other resummation techniques will 
give the same result, at least to leading order of the pinch-singular terms 0. 

Our new strategy is to suggest that in a general kinematical region where pinch singularities 
occur we have to solve simultaneously both KB equations for the propagator. It differs from 
the usual approach, where one treats the sum of the SD equations, for small K, as a constraint, 
which effectively puts Q on the mass shell. But how can one use two different solutions for 
the same propagator? As will be shown in Section ^, at one loop level these contributions 
are relevant to different kinematical regions: in Boltzmann equations we must have fco < 2go, 
in the other case (called normal cut regime in the following) > 2go must be true. At one 
loop the exactly known complete solution is thus the sum of the respective KB equations. 
We propose to maintain this additivity feature also for the resummed propagator. We expect 
that the kinematical complementarity of the two solutions is a good approximation while the 
two-particle threshold is well separated from the regime where Landau- damping dominates. 

While in our approach no constraint equations play any role, in the leading order pinch- 
singular terms we can solve the SD equations with an ansatz implying the missing constraints 



(see Sections TA. and |4.2| ). In this way the usual solutions for the Boltzmann equations survive, 
we just have supplementary contributions to the propagator from the complete solution. 

Computing expectation values of composite operators involves integration over Q, which is 
the momentum associated with the relative coordinate. After this integration, more insight is 
offered into the physics of the different contributions. At one loop order the discontinuity of the 
Boltzmann term ~ 1/ (QK) gives a non-zero result for < - this is the Landau damping 
regime. The discontinuity of the other contribution ~ 1/(Q^ + — m?) is non-zero when 

> Am? where 2m is the threshold value for particle production - this is a normal cut regime. 
The solution of the Boltzmann equation in the relaxation time approximation generates a finite 
imaginary part to k^] this is why a quantity, conserved at tree level, can decay exponentially 
due to radiative effects. In a kind of similar relaxation time approximation, higher loop pinch 
singularities provide imaginary and real contributions to the denominator Q'^+K'^ / A—m? , which 
results in a shift of the threshold location and in a "smoothing" of the threshold singularity. 

In this paper, the effects described above in qualitative terms will be demonstrated in the 
cubic scalar model. Beyond its simplicity it shows some similarity with fermion-gauge systems 
as far as the perturbation theory is concerned. We choose a simple observable, {[ip'^{x), ^9^(0)]); 
this is the simplified analogue of ([j^(x), j'^(O)]), which plays an important role in the theory of 
dilepton production |]TU|. In Section ^ we introduce the model and our observable. In Section ^ 
we perform a thorough one-loop calculation, identify the different kinematical regions, and 



2 



examine the way SD equations and KB equations reproduce the one-loop result. In Section | 
we turn to higher loops. We show the origin of pinch singularities, and apply KB equations 
to resum these singularities in different kinematic regions. We also solve the corresponding 
equations in relaxation time approximation. In Section |^ we make some comments on the 
relevance of our results to the dilepton production in fermion-gauge systems. Section |^ then 
closes the paper by presenting our conclusions and an outlook. 



2 The model 

The cubic scalar model has the Lagrangian 

C = \id^r-lm^^'-^^^ (3) 
We are interested in the equilibrium expectation value of the composite spectral function: 



(4) 

lPI,{ip)=0 



This quantity corresponds to the expectation value {[j^i{x),j^{0)]) in the fermion-gauge sys- 
tems, which plays a central role in dilepton production |T^. The subscripts IPI and {(f) = 



will be omitted in the following. The quantity of phenomenological interest is the discontinuity 
of iR{K) , which is the Fourier transform of 

zR{x) = Q{xo){[v\x),^\0)]). (5) 

R{x), on the other hand, is the linear response function of {ip'^{x)) to the mass modification of 
the Lagrangian characterized by the linear operator V: 

Cr = C-^V^' 6 {^') (K) = ^R{K)V{K) + . . . , (6) 

where Sf = f — /|-p=o for any quantity. We will define {ip'^{x)) via point splitting; in the 
Keldysh and in the R/A formalisms (cf. Appendix ^ we should compute 

[f'^ix)) = ^\im[{(f{x)(f{y)) + {if{y)ip{x))] = ^lim[iG2i{x,y) + iGi2ix,y)] = limiGrr{x,y). 

(7) 

In Fourier space 

y,2^ (K) = Jclxe'^'^lim Jd;pd:qe-'P''+'''HG„{p,q) = J dQ iGrr{Q, K), (8) 

{G{Q,K) denotes the Wigner-Fourier transform of G{x,y), which has the usual two- variable 
Fourier-transform G{p,q), cf. Appendix ^). We introduce the quantity R{Q,K) defined 
through the relation 

5iGrr{Q,K) = ^R{Q,K)V{K) ^ R{K) = J dQ R{Q , K) . (9) 
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3 One-loop perturbation theory 



In leading order the only interaction we have to take into account is In the Keldysh 

and R/A formalism (cf. Appendix 0) this piece of the Lagrangian has the form 



--V [^l - 



Linear response to V yields 
5G„{Q,K) 



(10) 



V{K), (11) 



where means free equilibrium propagators. The coefficient of V identifies —iR{Q,K)/2. 
Its discontinuity reads 



BisciR{Q,K) = 2 

ko 



q\Q + ^)^Gl{Q - y) - q\Q - ^YGl{Q + |) 



Inserting the free equilibrium expressions (cj^ = + m? 



^°(Q) = 27rsgn(go)5(go' - ^q), ^C(Q) = 2vrsgn(go) (^^ + ^(^o) j 5{ql - uD 
into this equation yields 



DisczR{Q,K) = 87r'sgn{qo+'^)sgn{qo-^){n.-n+)6{{Q+^?-m')6((Q-^?-m'' 
ko 2 Z 2 



(12) 



(13) 



(14) 



where n± = n(go ± fco/2). 

We first analyse the product of the delta functions (for details, see Appendix One can 
write them in various forms: 



D := 5iiQ + |)2 - m')Si{Q - - m') = 5{2QK)5{Q' + ^ - m') 
= Wi '^^go - Y'Sn(A;o)2)5(g - ^2), 



(15) 



where Z = \J {K"^ — Aw?) / {k^ — k'^x'^), which is real when > Aw? or < (see Ap- 
pendix 0) . 

The last form carries an important message concerning the relative magnitude of go and k^. 
For K"^ < it requires that |go| > |fco|/2, while for K'^ > 4m^ we find |go| < k/2 < |fco|/2 (see 
Appendix ^). It has the consequence that in the present circumstances one can write: 



sgn(go ± y ) = QiK' - 4m2)(±)sgn(fco) + e(-K2)sgn(go), 



(16) 



which implies, when inserting back into ([14|) , that 

47r2 



DisciR{Q,K) 

ko 



r-. ^{(lo - ^sgn{ko)Z)S{q - ^-^Z) x 

^{k^ - k^x^){K^ -4w^) 2 2 



- 4^2) (l + - go) + n(| + go)) + &{-K') Lq, - |) - n(go + |; 



:i7) 



The Q integration is performed using this form (see Appendix P), and one finds 



kZi/2 

/^QDisczi?(Q,i^) = ^^[ 0(^^-4m2) / dy U + n(}-^ - y) + n(}^ + y)\ + 

-kZi/2 ^ ^ 

(18) 



2Q{-K') J dy U(y-^)-n(y + ^) 



kZi/2 



which agrees with the results of earher calculations |Tl[ . 

The limiting case usually investigated is that the spatial variation of the external source is 
small, i.e. k <^ m,T. Then, in the first term of eq. (|TB]), the variation of y is of order k, thus 
y ^ |A;o| and we can neglect the y dependence of the integrand. In the second term k > |fco|; 
so the small-fc approximation yields Zi ^ 2m/\K\ S> 1. Then y ^ k > ko and we can perform 
an expansion in power series with respect to ko. Finally we obtain 



/cfQ I)isciR{Q,K) ^ 

J ko 4:TT 




^ ko km. 



■ (19) 



3.1 The Landau damping region 

The approximate expression (|l^) can be derived without the detailed computation outlined 
above. If < and k <^ m, T, the plasma excitations behave as on-mass-shall particles 
because = m'^—K'^/4 ^ m^. At high temperatures all components of Q are of the order of the 
temperature. In this regime it is a good approximation to neglect K in (0) and ( pTS]) whenever 
it is possible, and to perform an expansion in power series with respect to the components of 
K when it is necessary. This yields 

Disc«i?(g, K) ^ -Sn^ko ^^^^6{2QK)6{Q^ - m') = -87r^ko^6{2QK)6{Q^ - m^), (20) 
ko ago du 

since dn/dqo is even in qo- After the Q integration one readily reproduces the Landau damping 
piece of (0). 

Although this is a straightforward consequence of our approximations, this result is still a 
bit surprising. The assumption we made is quite plausible: the variation of the background 
field K is small as compared to the typical scale of the loop momentum, which is T at high 
temperatures. This was implemented by requiring Q ^ K for all components. This assumption, 
however, just reproduces the Landau cut region, while the normal cut is missing. 



3.2 The normal cut region 

Let us examine where we have lost the normal cut, and what are the appropriate assumptions 
that lead to its approximate form in (p!9|), in a similarly direct way as the Landau cut is 
reproduced. 

The answer is easy to find when we recall the earlier statement of the exact analysis that in 
the > 4m^ regime |go| < k/2 < \ko\/2: that means that the assumption Q ^ K is not valid 
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for all components. In the normal cut regime the time-like component of the loop momentum 
is small, even smaller than the spatial variation of the background field. 

This observation suggests a complementary approximation scheme: we neglect the (spatial) 
k dependence, not assuming, however, that ko < k. In the product of the delta functions (|T5|) 
we can then write 

D = 6{2QK)6{Q' + ^-m')-^ 6{2qok,)5{Q^ + ^-m^) = _l_5(g^)5(M _ ^2 _ ^2y (21) 

Therefore, in the leading order, we can assume go ~ as well. We also see that we immediately 
obtain k^ > 4m^ from the second delta function. The approximation for the discontinuity of R 
then reads 

471-2 / u \ 1.2 

Disczi?(Q, K)^ — [l + 2nA Siqo)SQ - - q'). (22) 

This expression reproduces, after Q integration, the normal cut part of the approximated result 

m- 

It is interesting to observe that the approximations leading to Landau and normal cuts 
concern different factors of D in (|l^): using 6{Q'^ + — m?) for small K yields that Q 

is on the mass shell, while using 6{2QK) for small k yields go ~ 0. While both assumptions 
seem to be well justified, neither of them reproduces the complete discontinuity of R, just 
complementary parts of them. To understand this phenomenon better, we now turn to the 
Schwinger-Dyson equations. 

3.3 Schwinger— Dyson equations 

The general form of SD equations in Keldysh and R/A formalism can be found in Appendix |^. 
At one loop perturbation theory we need the self-energies coming from the interaction La- 
grangian ([T0|) : 

UariQ, K) = n,,(Q, K) = PiK), Uaa = 0, (23) 

which has to be inserted into (|8TD of Appendix 0. For a linear response we insert free propa- 
gators into its right-hand side. We are interested in the rr component (cf. (^) 

(g ± |)2 - m^) GrriQ, K) = V{K)Gl{Q t f )• (24) 

The two equations, like the general ones, are related by complex conjugation and K —K 
transformation. 

Taking the sum and the difference of the two equations we find 



2QKGrr{Q,K) =V{K) 



Gl{Q-^)-GliQ + ^) 



(26) 



This yields two solutions for Grr{Q,K). Since the two SD equations are not independent, the 
two solutions should imply the same physics. To see this we recall that (5°^(Q) ~ S{Q^ — m^). 
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therefore 



g2 + ^ - m2 



2QK 



(27) 



On the other hand the equivalence is not complete, since we lose something from the analytical 
structure. If we treat these equations as differential equations with initial condition (i.e. we 
apply the Landau prescription), and we take the imaginary part, the two equations give different 
results. We use the following form of the free rr propagator (o;^ = + m^): 



Disc iGrriQ, K) = 47r2p(ir)sgn(fco) (1 + n{uj+) + n{cu.)) S{Q^ + -r - m^)5{2QK), 
T)isciGrr{Q,K) = 47r^P(ii')sgn(go) (n(cj_) - n{Lj+))6{Q^ + -r ' m^)K'^QK), (29) 



respectively, where ijJ± = Ci;(q±k/2). Here the first form corresponds to the normal cut regime, 
the second to the Landau-cut regime. Although formally there seems to be nothing that could 
prohibit applying, for example, the second equation to the K"^ > 4m^ region, the result will be 
incorrect. 

If we compare these formulae to the ones coming from direct perturbation theory (see 
(p!^ and ([T5|)) we can see that the difference is that that in the exact formula there stands 
sgn(go ± fco/2), which is substituted by sgn(fco) and sgn(go) in the first and the second form of 
the above equations, respectively. The reason for this is simple: when we take for example the 
difference of the two equations, disappears, and so the corresponding ikoe will be missing 
in the later calculations, too. Therefore we implicitly assume that the 2iqQe term (coming from 
2QK) dominates it: i.e. go > This implicit assumption prevents the final result to be 

applicable everywhere. 

Usually one considers this pair of equations as a differential equation ((|26|), the Boltzmann 
equation) and, assuming K <ti Q, a. constraint (p5|). But this strategy, as we have seen, misses 
the normal cut part from the very beginning. From our analysis a slightly different picture 
seems to emerge: here both equations are differential equations to be solved, and the complete 
solution is the sum of the two solutions. We propose to apply this strategy also for the non- 
perturbative, non-linear response regime: we should solve not just the equation coming from the 
difference of the pair of SD equations, i.e. the Boltzmann equation. We have to treat the sum 
of the two equations as a differential equation, since this is the approach that makes explicit 
where the potentially important (especially in the zero-mass case) normal cut contribution 
comes from. The complete solution then can be composed as the sum of these partial solutions. 

But where are the constraint equations? In our case the right-hand side, and so the solution 
was explicit, there was no need to add another equation. At higher loops or in the non- 
linear response regime, as we will see in the next section, concentrating on the IR sensitive 
pinch-singular region, we can still have a consistent ansatz that satisfies the (by now missing) 
constraint equation. 





and write for the solution of (p5|) and ( pG]) 
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As a last remark in this section we emphasize the importance of the cut contributions in free 
systems. Let us consider the evolution of / (i^x(y9^)(x, /cq) = (v^^)(k = 0, /cq) in a free system, 
without any background field, but with specified initial conditions. The relevant equations are 
the following (cf. (|2|), (|2|)): 

2qokoGrr{Q, ko) = 0, {q^+'j- GrriQ, K) = 0. (30) 

Taking into account only the first equation we would conclude that Grr{Q,t) is constant, and 
thus (v3^)(k = 0, t) is constant, too. This would mean that (</3^)(k = 0,t) is a conserved quantity 
- which it is not, in reality. Analysing the second equation (cf. |Tl|, 0), we find that it yields a 
power law time dependence, that is ((y9^)(k = 0,t) goes to zero, with some power of time. This 
again shows that the Boltzmann equation, alone, cannot provide the complete treatment even 
in the simplest case. 



4 Higher loops and pinch singularities 

Higher loops are, of course, suppressed by powers of the coupling constant; therefore they 
become important only if the one-loop result is zero or there are IR divergences. The imaginary 
part of the one-loop result, as we have seen, is zero if < < Arn^; at two loops, however, 
this restriction is relaxed and non-zero imaginary part is obtained for all K. As far as the IR 
divergences are concerned, in our case the pinch singularities [0, |^ are the primary sources of 
them. 

If the product of a retarded and an advanced propagator appears in some objects with 
the same momentum, then the poles of the propagators pinch the real axis as e — > and 
yield a non-regularizable double pole. If the two momenta are not equal, we expect a singular 
behaviour of the object as the difference of the momenta vanishes. Using the fact that GrGr 
is not singular and that i{GR — Ga) = Q, we can write 

GniQ + y )G'A(g - y ) = GniQ + y) i0iQ - y ) - non-singular. (31) 

li K ^ the free spectral function g{Q) ~ 5(Q^ — itl'^) puts Gr on-shell: this is the core effect 
of the pinch singularities. On the other hand, Grr ~ Q, thus similar singularities appear when 
we multiply Gr^a and Grr- We have to treat this case in complete analogy with the GrGa 
product. 

Pinch singularities appear in all higher order diagrams containing two-particle intermediate 
states 0; to see the net effect we have to sum up all the diagrams of this type. In the hnear 
response case, these diagrams form ladders. 

Which are the kinematic ranges where we have to count with pinch singularities? From the 
above generic form 

GniQ + y ) wiQ - y ) - " 2 Q^ + KyA-m^' ^^^^ 
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There are therefore two sensitive regimes: the vanishing 2QK or the vanishing Q"^ + — iv? 
regime. 

To resum the pinch singular diagrams we can use any vahd equation where the left-hand 
side is proportional to 2QK or + — m^, respectively, and the right-hand side is non- 
vanishing in the 2QK = or + — m? = limit 0; the difference in the exact choice 
of the equation just gives a subleading effect. In linear response theory we need a linearized 
right-hand side. Exactly of this form are the linearized SD equations (|8TD of Appendix 0. 
Starting from it, and linearizing around equilibrium in the form 

G{Q, K) = G'''{Q)S{K) + g{Q, K), n(Q, K) = ff«(g)^(ir) + 7r(g, K), (33) 

one finds the following equations for the deviations: 

2QKgrr = Haa {p^r " G^a) + ^tadar - ^aa9ra + 



ra I 9rr 



+ ^ - 2gra = {G;, + G+ ) + (n+ + H" ) g^ 



2QKgra = -Kar [G;, " G'+ ) + (n+ - g^a, (34) 



where we have omitted the (Q, K) arguments of the functions and we have introduced the 
notations = F^'^{Q ± K/2). The V dependence is hidden in TT^a and iiar 

T^ra V + -Kra, T^ar^V + Hav (35) 

The self-energies (apart from the term containing V) come from loop diagrams, i.e. they 
consist of propagator products supplemented with some integration. In a linear response, all 
propagators but one are in equilibrium, so they are functions of space-time coordinate difference. 
Therefore the general structure reads in Fourier space 

7r(p,g)~ J dkdilC{p,~q,k,-i)g{kJ), (36) 

where we divided the momenta into incoming {p,k) and outgoing {q,i) ones, but this is just 
a convention in the signs. Since JC contains equilibrium propagators only, the momentum is 
conserved, and we can write 

7r(p,g)~ Jdi)C'{p,-q,i)g{p-q + i,i). (37) 

The Wigner transformed function reads 

tt{Q,K)^ J dilC'{Q,K,£)g{K + e,K) = J d£ lC"{Q, KJ)g{£, K). (38) 

In general, we have to sum over the R/A indexes, too. From this form we see that the linearized 
SD equations are diagonal in K. 
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4.1 The Boltzmann equations 

Instead of using the full set of SD equations in (0), we can make some approximations in order 
to construct equations that can be treated more easily. According to our earlier analysis, we 
will concentrate on the summation of pinch singularities, which covers two regions: the small 
QK region and the small + K'^/A — m? region. The approximation procedure will be to 
power-expand everything with respect to the powers of these small parameters. 

We start with the Boltzmann region: here go > ^o/2 and the differential equations are the 
second and fourth of (|3^). First we examine the equations for g^a'- 



2QKg,a = liar {c^a " G^a) + (n+ " H" ) g,,. (39) 
If go > ^o/2 and in the small QK case the following approximation is valid 

^™ ^ + K^/A-m''±QK + iqoe ^ + K^A - + iq^e ^ ^^^^ 

therefore the difference of these functions is zero up to 0{QK). Then the equation for g^a is 
homogeneous (no mixing with grr)-, therefore a valid ansatz is to assume that g^a = 0. In a 
similar way we can assume that gar = 0. This is a generally accepted approximation: the usual 
Boltzmann equation never leads to non-zero gra and gar- 
Having said this, the equation for grr reads 

2QKgrr = TTaa (C^, - ) + T^arG^ - 71^6'+ + (^11+ - II^^) g„. (41) 

We approximate the free Grr for small 2QK and go > A;o/2 as (here n± = n(go ± /co/S)): 
tG% = 27rsgn(go±^) + n±) 6{Q' + ^-m'±QK) ^ 27rsgn(go) + n±) 6{Q' + ^ 



m y. 



while the difference between Gra and Gar, as we have seen in (|40|), reads as 



(42) 

;en in (|40|), reads as 

G-, - G+a = 2msgniqo)6iQ' + ^ - m') + 0{QK). (43) 

Therefore all explicit terms in the equation are proportional to sgn(go)(5((5^ + -ft'^/4 — m^); a 
consistent solution, therefore, comes by writing the ansatz 

igrriQ, K) = 27rsgn(go)5(g2 + ~ ~ m^)5n{Q, K). (44) 

Therefore, without making use of the other equation arising from the sum of the two original 
SD equations, we are led to a form consistent with this constraint. This is now, however, 
somewhat weaker statement, it is just an ansatz, and there may exist other solutions that do 
not respect this ansatz. 

If we are interested in the small k limit of the result, we can make other simplifications as 
well. Since in this region < 0, it also means a small ko, and we can thus neglect on the 
right-hand side. We find finally 

2QK5n = -Kaa + (vTar - VT,,) ^ + n(go)) + (n°, - n°J5n. (45) 
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This is the hnearized Boltzmann equation when we expand vr's in terms of 6n. 
In our concrete model (Rf), using the ansatz, we find 



Tin 



-th^ J dL q{L)q{Q - ^) [(^ + n{q^ - ^o)) 6n{^o) + + n{^^)^ 6n{qo - io) 

TTar - TTra = -ik^ J q{L)q{Q - L)[5n(^o) + ^^^(go " ^o)]- (46) 

The equihbrium value of the self energy reads 

Kr - Ka = -^h' J g{L)g{Q - L)(l + n{io) + n{qo - 4)), (47) 
and thus for the integrand on the right hand side we find 

6n{io)in{qo) - n{qo - 4)) + 6n{qo - io){n{qo) - n(4)) + 5n(go)(l + n{io) + n{qo - io)), (48) 
while the complete equation is 

2QK6n = -th^ J dL g{L)g{Q - L) 6[ (1 + n{io)){l + n{qo - io))n{qo) 

-n{£o)n{qo-eo){l + n{qo))]. (49) 

This is, indeed, a linearized Boltzmann equation. 

This expression is actually zero, since = m?, = m? and [Q — Vf' = m? cannot be true 
at the same time. In our model, therefore, there is no one-loop correction to the Boltzmann 
equation. To get a non-zero right-hand side, we can consider a model with different masses 
where nii > 1712 + rris is true, or we go to higher loops in the present model. 

4.2 The normal cut region 

Let us now repeat the same analysis for the cut region, where ko > 2qo and the first and the 
third differential equations of (0) are relevant. First, we take the equation for gra 



{q^ + ^- 2(7.a = TTar (G" + G+ ) + (n+ + B" ) 



(50) 



In this region the following approximation is valid 



1 ±1 K'^ 

" Q^ + Ky4-m''±QK±tkoe ^ QK + ik^e + + — -^^)^ (51) 

and so the sum G^^ + G^^ ~ 0. The equation is approximately homogeneous in g^a] we can 
therefore assume, as before, that gra = Qar = 0. 
Then the equation for g„ reads 

( + ^ - J 2grr = TTaa (G'^ + ) + TlarG'^ + TlraG^ + (n+ + B^^) grv (52) 
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We now approximate G„ for small + K^/A — m? and for kQ/2 > a.s 

^ 27rsgn(A;o) f ^ + ± go) ) SiQK) (53) 



2 ' 2 



and, according to the previous approximation: 



Gar + Gta ^ + ^ ^ ^ = -2'KisgTi{ko)5{QK) ■ (54) 

So finally all known terms are explicitly proportional to sg'a{kQ)5{QK), and we can have, 
similarly to the Boltzmann case, an ansatz that now reads as 

iQrriQ, K) = 2nsgn{ko)S{QK)f{Q, K). (55) 

This agrees with the solution of the difference of the SD equations, if we have considered it as 
a constraint; but now it is just an ansatz. 

In the small k region 5{QK) f» l/\ko\8{qQ), we can therefore substitute go ~ 0. Then the 
only non-zero momentum components are /cq and q, and we collect them in a single four-vector 
S = (/co/2,q). Using spatial rotation symmetry, liraiQ — K/2) = Ilar{S), and we can write 

(S^ - m^)2f = TXaa + (vTa. + T^ra) + n{so)^ + 2Iiar{S) f . (56) 

This is a linear but non-local wave equation for /. On the right-hand side we have to use the 
fact that Q"^ + — is small, i.e. S"^ ^ w? . 

In our concrete model, at one-loop order 7raa(/co/2, q) = IvciHariS) = at 5*^ = w?. The 
second equality was proved before. In order to prove the first one, we notice that, according to 
the ansatz: 

vraa(Q, K) = -2ih^ J dL g{Q - ^) + ^(^o - 4)) sgn{ko)27T 6 {LK)f{L, K). (57) 

When k = go = then 6{2LK) ~ and thus g{Q — L) would force Q — L on the mass 

shell with go — ^o = 0; which is impossible. Thus, at the one-loop order, what remains in ( |56D 
(writing the V dependence explicitly) is 



/I \ h"^ r li^f 1 /I \ 

(S^ - 'n')f = <(S)f - (5 + n(^o)) - + (5 + »(.o)) P(2.„). 

(58) 

We see that the treatment of the two regimes (the Boltzmann and the normal cut regions) is 
quite symmetric; the only difference comes from whether we consider the sum or the difference 
of the SD equations as the relevant dynamical differential equation. In the IR sensitive regions 
we have for the solution an ansatz that satisfies the other equation as a constraint. 
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4.3 Relaxation time approximation 

Relaxation time approximation for Boltzmann equations is well-known. It means that we 
approximate the right-hand side of ( ^9]) as —2iqQT B^n. If this approximation is valid, then 
5n{t) ~ exp(— rst). In our illustrative model, as we have seen, Vb = 0{li^), since it is zero at 
the one-loop level. 

This approximation scheme can be extended also to the normal cut regime. We assume 
that the operators acting on the right-hand side of the relevant eq. ( [58| ) can be represented 
as a diagonal linear operator. After rearrangement on the left-hand side, modified parameters 
(damping rate and mass) will appear. Then reads as 



-r '-J- c \ 2 2 
-q -m^ff 



f=y- + n{^)y{K). (59) 



At one- loop order the imaginary part = , just like in the Boltzmann case; but, in general, 
it is non-zero at higher loops. Its solution is almost the same as the one- loop level result, just 
the finite imaginary part requires some care. The ansatz (^) yields 

= ^*(*' (5 + "(f)) (.„ + ,r,)V4-q^ - 

When compared with the definition ig„.{Q,K) = R{Q, K)V{K)/2, we find for / dQiR, after 
evaluating the go integral. 



1/1 k \ 7 



After performing the integration and dropping the infinite piece (mass renormalization), we 
finally arrive at 



The discontinuity, calculated as twice the imaginary part, reads 



<tQ^RiQ,K) = ^ 1 + 2n(:^) f ^° ^ -mj,,. (62) 



1 / ,koA 2 (k^ . koTA , ^ 



where = '"^e// + and 

1/2 



F(x, y) = ^[x^ 4^^^^) . (64) 



Let us first analyse some properties of this result 
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If Fc = then, because of y - 
(p!9|) with the substitution m — > M . 



0) = 0(x)a/x, we get back the normal cut part of 



Using perturbation theory without resummation, we see IR divergences (coming from 
higher loops) if m = 0. This is because ~ /i^, thus we can compute power series in 
M^, which yields 

—j^ + Oih"). (65) 



2F 
ko 



pert 



For = 



2F 



Af2=o 



yields a divergent (although integrable) behaviour. 
At A;o = 0, M ^ we get 



kn 



1/2 



A:o-*0 



(66) 



2F 

ko 



ko-*0 



AM 



sgn(fco 



(67) 



This is a finite (and non- analytic) result. 



While the free result exhibited a non-analytic square-root behaviour near the threshold, 
higher-loop corrections smoothen out this non-analyticity. A non-analytic behaviour im- 



plies power-law real time dependence [ITT] , |T2[, which now acquires an exponential damping 
multiplicative correction, just like the solution of the Boltzmann equation. On the other 
hand, at ko = 0, there appears a new non-analytic contribution. 

In order to see these effects we plotted the function 2F/ko in Fig. [l|. Curve a shows the 



1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 





a 



b /' 



/ c 



1 I I I L 



Figure 1: The function F/ko (cf. 



1, Fe = 0, c: M = 1.4, Fe = 0.1, d: M = 0, F 



for different parameters: a: M 
0.1. 



0, b: M 



M = Fc = result; b is the free result with m = 1; c shows how it changes at higher 
loops: it acquires a finite shift in the mass as well as a damping that smears out the threshold 
non-analyticity (in the plot we have used M = 1.4, F^ = 0.1). Finally curve d shows the 
M = 0, F, ^ case. 
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5 Comments on the fermionic case 



Let us make some comments on the fermionic case. In our simple model, in relaxation time 
approximation, the only change in the final result, compared to the free case, was a modi- 
fication just in the threshold function Q{x)^/x F{x,y), cf. (|63|) . This seems to be the 
consequence of the relaxation time approximation, not depending on the specific model, so we 
may assume that this feature generalizes to other theories, too. In case of dilepton produc- 
tion we have a fermionic system with approximately zero masses. The rate of the dilepton 
production is proportional to {[j^, j'^]){K) |T^; when computed perturbatively, this is very 
similar to our observable ([v?^, ip'^]){K) = J dQ Discfcg iR{Q, K). At k = the production rate 
is (b/x, j^])(^o)/^o ~ tanh(/5/co/4); if we assume the same behaviour for the higher loop pinch 
singularities in the fermionic theory as they appear in our model, we obtain changes in the 
production rate qualitatively similar to Fig. B Here curve a shows the free result, curve h 




Figure 2: Expected result of ([j^, j'^])(/co)/^o- Curve a is the one loop result, curve b is the 
result after resummation. Curve c indicates the size of the 2-loop corrections. 



shows the expected result of the resummation, applying the present analysis. Curve c indi- 
cates the size of the 2-loop perturbative corrections. This figure suggests that, although the 
two-loop correction is divergent in the fco limit, the resummed result can be convergent, or 
even zero. With the present approximation the discontinuity goes linearly to zero for vanishing 
momentum, and its slope is ~ rc/(4TM). Our approximate result is quite reminiscent of the 
numerical (MEM method) findings of Karsch et al. 



6 Conclusion 

Let us summarize the main results of the paper: 

• Refinement of the conventional treatment of Kadanoff-Baym equations. Starting from 
the pair of non-equilibrium Schwinger-Dyson equations, applying Wigner transformation 
and taking their sum and difference, we arrive at the KB equations. The conventional 
interpretation is the following: for mild variations of the average coordinates, one of the 
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equations is understood as a conventional first-order partial differential equation: the 
quantum Boltzmann equation. In the other equation one neglects the dependence on the 
average coordinates, as a consequence of which it becomes a constraint. 

We propose a somewhat different approach. Based on the thorough analysis of the one- 
loop contribution, and resumming the IR (pinch) singularities arising from higher loops, 
we argue that both equations have to be interpreted as dynamical differential equations, 
each one beeing relevant in complementary kinematical regimes. The complete solution 
is the sum of the two equations. Concentrating on the pinch-singular contributions, we 
suggested for both solutions an ansatz that satisfies the constraint equation of the con- 
ventional interpretation. In this way one of the equations is still a (quantum) Boltzmann 
equation; the other one, however, represents a new contribution, reproducing in a re- 
summed form the particle production cut. 

• Features of the new contribution. The new contribution is responsible at tree level for 
the appearance of the normal cut; it is IR (pinch) singular at higher loops. We have 
constructed the relevant equation near the IR-singular regime (see (0)). It is a wave 
equation for the quantity that is the analogue of the particle number variation in the 
Boltzmann case. 

• Solution in the relaxation time approximation. We can try to diagonalize the linear 
differential equations and keep the eigenmode with the smallest damping rate: this is 
the relaxation time approximation in effect. In the Boltzmann regime this procedure 
changes the infinitesimal ie term of the Landau prescription to the finite imaginary part 
iTs, associated with the presence of an effective term in the equation of motion, which 
is linear in the time derivative. In the normal cut regime there are two contributions: 
a real one, which is in fact a mass shift, and an imaginary one, which is analogous to 
the damping rate in the Boltzmann case. Computing the normal cut contribution with 
this approximation we can observe two effects: the shift of the threshold location and 
the smoothing of the threshold non-analyticity. If treated in finite orders of perturbation 
theory, both effects may lead to divergences. 

• Prospects for generalization. The present computation was performed for the case of the 
cubic scalar theory, and direct computations in other more relevant models are required 
before one can have quantitative predictions: this is a project for the future. Still, some 
findings of the present calculation seem to be robust enough to conjecture similar features 
to be valid in other cases as well. The "dual" interpretation of the KB equations, the wave 
equation form of the equation describing the cut contributions, the imaginary and real 
parts in the relaxation time approximation probably generalize to the fermionic case as 
well. Then we conjecture for the quantity {[jf_i,j^]){ko)/kl, relevant to the determination 
of the dilepton production rate, the behaviour indicated on Fig. ^ While two-loop order 
perturbation theory gives a divergent contribution (curve c) at threshold, the resummed 
result shows a much milder behaviour (curve b). This qualitative prediction displays 
remarkable similarity to recent numerical data [|I|. 



16 



Acknowledgements 



The author would hke to thank to D. Litim, A. Patkos and P. Petreczky for discussions and 
proposals. This work was partially supported by the Hungarian Science Fund (OTKA). 



A Appendix: The formalism 

In the Keldysh formalism for a complex scalar fieldQ we define the propagators 

zGuit, t') = e{t - 1') {m^Kt')) + 0(i' - 1) {^Kt')m) , ^G'i2(t, t') = {^Ht')m] 

where the expectation value is defined in general as 

{A) = TrpA, (69) 

where p is the density matrix. Formally we introduce two types of fields ($i and $2), which 
are representations of the original field on the different branches of the closed time contour. 
Their expectation values are defined as = iGij. 

These propagators are not independent, there is a relation between them 

Gil + G22 = Gi2 + 6-21. (70) 

In fact, there are only two independent expectation values, e.g. Gu and G21, since Gu and 
G22 are combinations of these two and the theta function. In equilibrium even these two are 
related through the KMS condition [0. 

While these propagators determine the free theory, in case of interactions we have to modify 
the interaction Lagrangian as (cf. |0]) 

Ci{<^)^Ci{^i)-Ci{^2). (71) 
In the R/A formalism we introduce new fields 

$1 + $9 

<fr = ^ , $a = $1 - $2. (72) 

Their time-ordered products (denoting ^i{x) by $j and $J(?/) by $'•) are expressible through 
the two-point functions as follows: 

Grr = ^ ($1$; + $1$^ + $2$; + $2$^) = 7(^11 + G12 + G2I + G22) = ^i^±^, 

At 4 2 

Gra = ^ ($1$; - $1$ 2 + '^2<f'l - $2$ 2) = ^(^11 " G12 + G21 " G22) = Gn - Gu, 

Gar = ^ ($ l$'l + $ 1$ 2 - '^'2<f'l " $2$ 2) = ^(^11 + G12 " G2I " G22) = Gn - ^21, 

Gaa = -% ($l$'l - $1$2 - $2$'l + $2$2) = ^11 - ^12 - G21 + G22 = 0. (73) 



^For fcrmionic fields the formalism is analogous, except for some appropriately placed minus signs, see 
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In a non-equilibrium situation the propagators depend on the two space-time coordinates 
separately, and one can define several integral representations for the propagators. The double 
Fourier transform of a function of two arguments is defined as 

fip,q) = ldxdye"yfix,y), f{x,y) = J dpdqe-'P^+"^yfix,y), (74) 

where dp = d^p/ (27r)^ in D dimensions. Note that we define the momentum of the y vari- 
able with a relative negative sign! We use the same symbols for the function and its Fourier 
transform, the type of the argument defines which function we are dealing with. The Wigner 
transform / of a function / with two arguments is defined as 

/>, X) = /(X + ^, X - |), fix, y) = /> - y, ^), (75) 

i.e. it understands the function to depend on the relative and average coordinates. Fourier 
transformation with respect to these coordinates yields 

/(g, K)^j dudX e^'3"+'^^/(i., X)^ j dxdy e^Q(-.)+^^M/2/(^, y) ^f^Q+^^Q_^^^ 

fip,q) = fi^^-^,p- q)- (76) 

The derivatives with respect to x and y transform as 

idJix, y) p/(p, g) ^-^^^ (Q + ^)/(Q, X), 

^^,f{x, y) -g/(p, ,) (-g + ^)/(g, K). (77) 

Later on we will use the Schwinger- Dyson equations for propagators. It is used to resum 
bubble series 

Gab = Gab + Gac^cdGdb + . . . = Gab + G'^Jlcd.Gdb = Gib + Gac^cdG^b- (78) 

The indices refer to space-time coordinates as well as to internal indices (e.g. Keldysh indices). 
To have a diff'erential form we apply the free time evolution kernel on both sides. In the R/A 
formalism these kernels satisfy "twin" equations, which differ only in the variable to which the 
derivatives are referring to: 

{dly + m')G%{x,y)^0 and {dly + m')Gl{x,y) ^ {dl^ + m')Gl{x,y) ^ -S{x - y). (79) 

Using these, we find the SD equations: 

-{dl + m^)Grr ^ jdz {UaaGar + ^arGrr): -{d^ + m^)Grr ^ j dz {Gra^aa + Grr^ra) 
-{dl + m^)Gra ^5{x-y) + j dz UarGra, -{Oy + m^)Gra ^ S{x - y) + j dz Gra^ar: 

-{dl + m^)Gar ^5{x-y) + J dz n,„G„„ -(d^^ + m^)Gar ^ 5{x - y) + J dz G^.n.^. (80) 
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Here we have used the fact that Gaa = is consistent with SD equations only if 11^^ = 0. After 
the Wigner-Fourier transformations the equations in ( [80| ) look a little more complicated: 



K\2 
2 ' 



m 



') GrriQ, K) = Jdk 



K\2 
2 / 



I<\2 
2 I 



{Q 

\Q-fr 



m") Gra{Q,K) 
m^) GraiQ,K) 
Gar{Q,K) 
Gar{Q,K) 



T^aaiQ +iK- k)Gar{Q -^.k) + 

+ TlariQ + I,K- k)G„{Q -^,k) 

TlaaiQ -l,K- k)Gra{Q + ^ , k) + 

+ UraiQ -iK- k)GrriQ + ^) 



^(K) + JcfkUariQ +iK- k)Gra{Q - 
S{K) + !dktiar{Q -iK- k)GraiQ + 
S{K) + !dk UraiQ +iK- k)Gar{Q ' k) 



K-k 
2 ' 

K-k 
2 ' 



k) 
k) 



S{K) + j<ik IlraiQ -l,K- k)Gar{Q + 



K-k 
2 ' 



^1) 



where S{p) = {2tt)^6{p) in D dimensions. The ±K equations are related to each other by 
complex conjugation and K — > —K substitution. 



B Appendix: Detailed computation of Section @ 

In the product of delta functions (|15D we observe that the first delta function + QK + 
irV4-m2) implies QK = m'^ - Q"^ - K"^ /A, which yields S{Q^ -QK + K^/A-m'^) 5{2QK). 
This, in turn, means that one can set QK = in its coefficient, and therefore we can write 
5{Q'^ + QK + K'^ /A — m?) — >• + K'^ /A — w?). We obtain in this way the second form of 
(plSf ). In order to get the third form the following steps are done: 



D 



S{2qokQ 
l^ol 



K^ 



2qkx)6{qi - + — - 



1 







2|A;o 



qkx 
ko 



2 |/co 



k'^x^ 



■5iq 







qkx. ^, n kl K"^ — Am? . 




^0 



A k'o- k^x"^ 
6{q-^fZ)+6{q+^^Z) 



qkx 1 

2 ,.2^21 ^(^0-—)- 



2 \k^-k'^x 
1 



^0 



2q 



2 2^ I /cq 



k^x^ 



-S{qo 



ySgn(/i;o)2)5(g- ^2). 



k'^x^ 



m 



12) 



In the third line we have introduced the quantity Z = \J{K'^ — Aim?)/{kQ — k'^x'^). In order to 
get the fourth line we used the fact that, since q is positive, only one of the two delta functions 
appearing on the third line can be satisfied, therefore q = \kQ\Z/2. The expression is non-zero 
if Z is real, i.e. when (fT^ — 4m^)/(A;Q — fc^x^) > 0. To examine this condition we split the K'^ 
range into two regions: 
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• K"^ > 4m^; then /cq > /c^ > k'^x^ (since \x\ < 1), so the numerator and the denominator 
are positive; 

• < 4m^; then the numerator is negative, the denominator thus has to be negative as 
well. This yields < k'^x'^, which implies 1 > |x| > |fco|/^ and fc^ > k^, i.e. < 0. 

Therefore Z is real if > Am? or < 0. 

We can examine the relative size of |go| and |A;o| from the last form of ([T5|). To find it we 
rewrite Z as 



2 K^-Arn^ + 

2^ = = \ i L 

2^2 Z^2 ,T-»2 Z^2,T<^2 



3) 



The numerator is always positive, thus 



• if is:^ > 4^2^ ^Yien \x\ < 1, 2 < 1, and one has |go| = < f < 

• if ^2 < then 1 > |x| > l^ol/A;, Z > 1, giving |go| = > ^ 
For the Q-integration appearing in 

47r2 



IM- 

2 ' 



Disc ii?(Q,K) 



r-, '^(go - ^sgn(A:o)^)5(g - x 

^(A;g-Px2)(ir2_4^2) 2 2 



-^0 



e(ir^ - Am') 1 + - go) + + go) + e(-ir2) n(go - ^) - n(go + ^ 



we write the integration measure as 



(84) 



dqQ dq dx q^. 



i5) 



The following change of variable is made: 



_ kx dy _ klk K'^ - 4^ 
'^'^y=Y ' di " ~Y\ {kl - k^x'f 



1 r [-kZ,, kZi], 



if K'^ > Am'' 



^ ^ 2 1 [-00, -kZi] U [A;Zi, oo], if K' < 



16) 



where Zi = Z{x = 1) = Jl — Am'/K'. Then we have 



87r3 



oo oo 

dqo jdq 

oo 



kZi/2 



Q{K'-A') / dy + Q{-K') / dy + Q{-K') dy 



kZi/2 



2g2 Hkl-k'x')^ 



(87) 



Putting the ingredients together we find the result quoted in (p!8|): 

kZi/2 



dQI)isciRiQ,K) = ^-^^ 
ko Snk 



e{K'~Am') J dy ll + n{\^-y)+n{\^ + y) 



-kZi/2 



2&{-K') J dy Uy-^-^)-n{y + \f 



(88) 



kZi/2 
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